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Abstract 

We consider the Cauchy problem for a n x n strictly hyperbolic system 
of balance laws 

ut + f{u)x = g{x,u), xeR,t>0 
u{0, .) = Uo e LI n BV (R; R"), 
* |Ai(ti)| > c> for aU i e {1, . . . ,n}, 
\\g{x,-)\\c2<Mix)ehl, 

each characteristic field being genuinely nonlinear or linearly degener- 
ate. Assuming that the norm of ||<7(a;, and ||iio||bv(m) are small 
enough, we prove the existence and uniqueness of global entropy solutions 
of bounded total variation extending the result in [1] to unbounded (in 
Loo) sources. Furthermore, we apply this result to the fluid flow in a pipe 
with discontinuous cross sectional area, showing existence and uniqueness 
of the underlying semigroup. 
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1 Introduction 

The recent literature offers several results on the properties of gas flows on 
networks. For instance, in [4j [5l [6l [8] the well posedness is established for 
the gas flow at a junction of n pipes with constant diameters. The equations 
governing the gas flow in a pipe with a smooth varying cross section a{x) are 
given by (see for instance [IT]): 
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{dp , dq_ _ a'(x) „ 
at ^ dx ^ "S^y 
dq , 9{^+P} ^ _a'(a;) 
dt dx a{x) p 

The well posedness of this system is covered in [I] where an attractive unified 
approach to the existence and uniqueness theory for quasilinear strictly hy- 
perbolic systems of balance laws is proposed. The case of discontinuous cross 
sections is considered in the literature inserting a junction with suitable cou- 
phng conditions at the junction, see for example [ll[5l[9]. One way to obtain 
coupling conditions at the point of discontinuity of the cross section a is to take 
the limit of a sequence of Lipschitz continuous cross sections converging to a 
in LI (for a different approach see for instance [7]). Unfortunately the results 
in [T] require Loo bounds on the source term and well posedness is proved on 
a domain depending on this Lcxd bound. Since in the previous equations the 
source term contains the derivative of the cross sectional area one cannot hope 
to take the limit a*^ a. Indeed when a is discontinuous, the Loo norm of (a^)' 
goes to infinity. Therefore the purpose of this paper is to establish the result in 
[1] without requiring the Lcxd bound. More precisely, we consider the Cauchy 
problem for the following n x n system of equations 

ut+ f{u)^ ^ g{x,u), xeR,t>0, (1) 

endowed with a (suitably small) initial data 

u(0,x) = ito(x), xeR. (2) 

belonging to LI n BV(M;R"), the space of integrable functions with bounded 
total variation (Tot.Var.) in the sense of [12]. Here u{t,x) £ M" is the vector of 
unknowns, / : £7 — > R" denotes the fiuxes, i.e. a smooth function defined on Q 
which is an open neighborhood of the origin in R" . The system ([T]) is supposed 
to be strictly hyperbolic, with each characteristic field either genuinely nonlinear 
or linearly degenerate in the sense of Lax [10| . Concerning the source term 
we assume that it satisfies the following Caratheodory-type conditions: 

(Pi) g : R X n ^ R" is measurable with respect to (w.r.t.) a;, for any u e il, 
and is w.r.t. u, for any x £ R; 

(P2) there exists a LI function M{x) such that \\g{x, ■)||c2 < M{x); 

(P3) there exists a function uj G L1(R) such that \\g{x, •)||ci < ^j^ix). 

Remark 1. Note that the LI norm of M{x) does not have to be small but 
only bounded differently from oj{x) whose norm has to be small (see Theorem 
[T] below). Furthermore condition (P2) replaces the Loo bound of the norm 



2 



of (7 in [I] ■ Finally observe that we do not require any Loo bound on lo 
other hand we will need the following observation: if we define 

£h = sup / uj{x + s) ds, 

by absolute continuity one has £/i — > as — > 0. 

Moreover, we assume that a non-resonance condition holds, that is the charac- 
teristic speeds of the system ([1]) are bounded away from zero: 

|A^(u)|>c>0, \f uefl, ie{l,...,n}. (4) 

The following theorem states the well posedness of ([1]) in the above defined 
setting. 

Theorem 1. Assume (Pi)-(P3) and If the norm of u; in L1{M.) is suffi- 
ciently small, there exist a constant L > 0, a closed domain!? of integrable func- 
tions with small total variation and a unique semigroup P : [0, +oo) x V ^ T) 
satisfying 

i) PqU — u, Pt+sU ~ Pt o PgU for all u,v €z V and t, s > 0; 

ii) \\PsU — PtvW^K^^) < l(^\s — t\ + \\u — for all u,v (zV and t, s > 0; 

Hi) for all £ V the function u{t, •) = PfUo is a weak entropy solution of the 
Cauchy problem (^-(^ and satisfies the integral estimates ^4^ , 

Conversely let u : [0, T] ^ D be Lipschitz continuous as a map with values 
in and assume that u{t,x) satisfies the integral conditions ^J^, 

Then u{t, •) coincides with a trajectory of the semigroup P. 

The proof of this theorem is postponed to sections [3] and [4l where existence 
and uniqueness are proved. Before these technical details, we state the applica- 
tion of the above result to gas flow in section [51 Here we apply Theorem [1] to 
establish the existence and uniqueness of the semigroup related to pipes with 
discontinuous cross sections. Furthermore, we show that our approach yields 
the same semigroup as the approach followed in [5] in the special case of two 
connected pipes. The technical details of section [2] can be found at the end of 
the paper in section [5l 

2 Application to gas dynamics 

Theorem [T] provides an existence and uniqueness result for pipes with Lipschitz 
continuous cross section where the equations governing the gas flow are given 
by 

dp , dq a'{x) 



. On the 
(3) 



Here, as usual, p denotes the mass density, q the Hnear momentum, e is the 
energy density, a is the area of the cross section of the pipe and p is the pressure 
which is related to the conserved quantities (p, g, e) by the equations of state. 
In most situations, when two pipes of different size have to be connected, the 
length I of the adaptor is small compared to the length of the pipes. Therefore 
it is convenient to model these connections as pipes with a jump in the cross 
sectional area. These discontinuous cross sections however do not fulfill the 
requirements of Theorem [TJ Nevertheless, we can use this Theorem to derive 
the existence of solutions to the discontinuous problem by a limit procedure. To 
this end, we approximate the discontinuous function 

by a sequence a/ g C°^"'^(R, M+) with the following properties 

{a^, X < -i^ 
Vi{x), a:e[-^,^] (7) 
a+, 2; > ^ 

where <fi is any smooth monotone function which connects the two strictly 
positive constants a~, a"*". One possible choice of the approximations ai as well 
as the discontinuous pipe with cross section a are shown in figure [T] 




Figure 1: Illustration of approximated and discontinuous cross-sectional area 

With the help of Theorem [1] and the techniques used in its proof, we are now 
able to derive the following Theorem (see also [7] for a similar result obtained 
with different methods). 

Theorem 2. // ||aj||£]^ = — is sufficiently small, the semigroups 
related with the smooth section ai converge to a unique semigroup P. 

The limit semigroup satisfies and is uniquely identified by the integral esti- 
mates [44\ ), 114 5\ ) with [/" substituted by (see Sectionl^ for the point £. = 0. 
More precisely let u : [0, T] — s- P be Lipschitz continuous as a map with values 
in £1(IR,M") and assume that u{t,x) satisfies the integral conditions i44^ , ( |4<?| i 
with fj" substituted by for the point ^ = 0. Then u{t, •) coincides with a 
trajectory of the semigroup P. 

Observe that the same Theorem holds for the 2x2 isentropic system (see 
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Section O 



at ~^ dx 




(8) 



a (x) 
a{x) p ' 



dt + dx 



In [3] 2x2 homogeneous conservation laws at a junction are considered for given 
admissible junction conditions. The situation of a junctions with only two pipes 
with different cross sections can be modeled by our limit procedure or as in [6] 
with a suitable junction condition. If we define the function ^! which describes 
the junction conditions as 



then it fulfills the determinant condition in [SI Proposition 2.2] since it satisfies 
Lemma [3l Here <I>(a, u) is the solution of the ordinary differential equation (|62p 
in Section O With these junction conditions one can show that the semigroup 
obtained in [6] satisfies the same integral estimate (see the following proposition) 
as our limit semigroup hence they coincide. 

Proposition 1. The semigroup defined in IBl with the junction condition given 
by (0) satisfies the integral estimates \44\I j with [/" substituted by C/" for 
the point ^ = 0. 

The proof is postponed to Section O 

Remark 2. Note that Proposition 1 justifies the coupling condition ^ as well 
as the condition used in |9] to study the Riemann problem for the gas flow 
through a nozzle. 

3 Existence of BV entropy solutions 

Throughout the next two sections, we follow the structure of [1]. We recall some 
definitions and notations in there, and also the results which do not depend on 
the Lcx) boundedness of the source term. We will prove only the results which 
in [ij do depend on the Loo bound using our weaker hypotheses. 

3.1 The non homogeneous Riemann-Solver 

Consider the stationary equations associated to ([T]) , namely the system of ordi- 
nary differential equations: 



(9) 



f{v{x))x = g{x,v{x)). 



(10) 



For any Xo G M, w G f2, consider the initial data 



v{xo) = V. 



(11) 
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As in [T], we introduce a suitable approximation of the solutions to (|10|) . pT|) . 
Thanks to (|4]), the map u t-^ f{u) is invertible inside some neighborhood of the 
origin; in this neighborhood, for small h > 0, we can define 



Throughout the paper we will use the Landau notation 0(1) to indicate any 
function whose absolute value remains uniformly bounded, the bound depending 
only on / and ||Af||Lr 

Lemma 1. The function ^hixo,u) defined in US\) satisfies the following uni- 
form (with respect to Xo (zR- and to u in a suitable neighborhood of the origin) 
estimates. 




(12) 



This map gives an approximation of the flow of (jlOp in the sense that 




(13) 



lim sup |$/,.(xo, u) — u\ = 0, 



lim \\Id- Du^h{xo,u)\\ = 



(14) 



Proof. The Lipschitz continuity of / 



^ and ([3]) imply 



<i>hixo,u)-u\ = |$,,(a:„,u)-/-i(/(7i))| < 0(1) 





Next we compute 



Du^h{xo,u) ^Df ^ f{u)+ / g{xo + s,u) ds 
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which together to the identity u = f ^ (/(u)) implies 



\\Du<p!,{xo,u) -Id\\^ \\Du<i>h{xo,u) - Df-^ {f{u)) II 

Df-^ f{u)+ / g{xo + s,u)ds 



< 



Df-'ifiu)) 



WDfiu) 



\\Dug{xo + s,u)\\ ds 



+ \\Dr'if{u))\\ ■ / \\Dugixo + s,u)\\ds 



< 0(1)£,, 0. 

Finally, denoting with Di the partial derivative with respect to the i component 
of the state vector and by ^h,e the £ component of the vector ^h, we derive 

DiDj^hA^o^u) = ^ ^DkDk'f^^ ^/(u) + g{xo + s,u) ds^ 
■ DJk{u) + / Digk{xo + s, u) ds 



Djfk' (u) + / Djgk' {xo + s, u) ds 



+ s, u) ds 



DjD Jkiu) + / D.jDigk{xo + s, u) ds 



so that 



D^^h{xo,u)\\ <0(1) il + j M{xo + s)ds \ <0(1)(i + ||M||li) < 0(1). 



□ 

For any Xo S M we consider the system ([T|) , endowed with a Ricmann initial 
datum: 



m(0, x) 



ui if X < Xo 

Ur if X > Xo- 



(15) 



If the two states ug, Ur are sufficiently close, let ^' be the unique entropic 
homogeneous Riemann solver given by the map 



= ^'(fi)(Uf) = '0n(CTn) O . . . O ■(/'l (cTl ) (uf ) , 
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where fi = (cti, . . . , (t„) denotes the (signed) wave strengths vector in R", [TU]. 
Here ipji j = 1, . • . , n is the shock-rarefaction curve of the j'^ family, parametrized 
as in [3] and related to the homogeneous system of conservation laws 

ut + I{u)^ = 0. (16) 

Observe that, due to H]), all the simple waves appearing in the solution of (fT6|) . 
(fT5|) propagate with non-zero speed. 

To take into account the effects of the source term, we consider a stationary 
discontinuity across the line x — Xo^ that is, a wave whose speed is equal to 0, 
the so called zero- wave . Now, given /i > 0, we say that the particular Riemann 
solution: 

u{t,x) = { ;f ^<^« vt>0 (17) 

is admissible if and only if Ur = ^h{xo,Ui), where is the map defined in 
(fTS]) . Roughly speaking, we require ui, Ur to be (approximately) connected by 
a solution of the stationary equations (fTO|) . 

Definition 1. Given h > suitably small, Xo € M, we say that u{t,x) is a 
/i-Riemann solver for H]), (|4]), (fT5|) . if the following conditions hold 

(a) there exist two states , which satisfy = ^h{xo,u^)', 

(b) on the set {t > 0, x < Xq}, u{t,x) coincides with the solution to the 
homogeneous Riemann Problem (|16p with initial values Ui, u~ and, on 
the set {i > 0, x > Xo}, with the solution to the homogeneous Riemann 
Problem with initial values w^, u^; 

(c) the Riemann Problem between and is solved only by waves with 
negative speed (i.e. of the families 1, . . . ,p); 

(d) the Riemann Problem between and Uj. is solved only by waves with 
positive speed (i.e. of the families p + 1, . . . , n). 

Lemma 2. Let G M and u, ui, U2 be three states in a suitable neighborhood 
of the origin. For h suitably small, one has 

\<^hixo,u)-u\ = 0(1)/ uj{xo + s)ds, (18) 

Jq 

\'i>h{Xo,U2)-<^h{Xo,Ui)-{u2-Ui)\ = 0(l)|?i2 -Ull / w(s)ds.(19) 

Lemma 3. For any M > there exist 5'i, h[ > 0, depending only on M 
and the homogeneous system \lb]] . sueh that the following holds. For all maps 
cjx^C^ (]R",R") satisfying 

||0||c2<M, \c^{u)-u\<h\, \\I~Dcj){u)\\<h'^ 
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w 














Ur 



Figure 2: Wave structure in an /i-Riemann solver. 

and for all Ui G -8(0, S[), Ur G B{(f){0), S[) there exist n + 1 states wq, . . . , Wn+i 
and n wave sizes cti, . . . , cr„, depending smoothly on ui, Ur, such that with pre- 
vious notations: 

i) Wo = ue, Wn+i = Ur; 

a) Wi = *i(cri)(wi_i), z = 1, . . . ,p; 

Hi) Wp+i = (t>{wp); 

iv) Wi+i ^ '^i{ai){wi), i^p+l,...,n. 

The next lemma establishes existence and uniqueness for the /i-Riemann 
solvers (see Figl2]). 

Lemma 4. There exist Si, hi > such that the following holds: for any Xq G 
h € [0, /ii], ui, Ur G i?(0, (5i), there exists a unique h-Riemann solver in the 
sense of Definition [IJ 

Proof. By Lemma[T]if /ii > is chosen sufficiently small then for any h G [0, hi], 
a^o G M the map u i— > <^h{xo,u) meets the hypotheses of Lemma [3l Finally 
taking hi eventually smaller we can obtain that there exists (5i > such that 
B{Q,6i) C B{0,S[)nB{^h{xo,0),S'i), for any h G [0,hi]. □ 

In the sequel, E stands for the implicit function given by Lemmas [3] and 01 

Q=E{h,ui,Ur;Xo), 

which plays the role of a wave-size vector. We recall that, by Lemma |3l is a 
C'^ function with respect to the variables ui, Ur and its norm is bounded by 
a constant independent of h and Xo- 

In contrast with the homogeneous case, the wave-size 6 in the /i-Riemann 
solver is not equivalent to the jump size \u^ ~ Ur\\ an additional term appears 
coming from the "Dirac source term" (see the special case U£ = Ur). 
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Lemma 5. Let Si, hi be the constants in Lemma^ For ui,Ur G B{0,Si), 
h€z[0,hi], set Q = E[h,U£,Ur'TXo\- Then it holds: 

\ui~Ur\ = 0(1) I |6| + / Lj{Xo + s)ds], 



|6| = 0(1) I \ni — + / uj(xo + s) ds . 



3.2 Existence of a Lipschitz semigroup of BV entropy so- 
lutions 

Note that as shown in [T| we can identify the sizes of the zero waves with the 
quantity 

a = I uj{jh + s)ds. (21) 

With this definition all the Glimni interaction estimates continue to hold with 
constants that depend only on / and on therefore all the wave front 

tracking algorithm can be carried out obtaining the existence of e, /i- approximate 
solutions as defined below. 

Definition 2. Given e,h > 0, we say that a continuous map 

u"-'' : [0,+oo) -^Lhoc (R,M") 
is an e, /i-approximate solution of (H])-® if the following holds: 

- As a function of two variables, it*^''' is piecewise constant with discon- 
tinuities occurring along finitely many straight lines in the x, t plane. 
Only finitely many wave-front interactions occur, each involving exactly 
two wave-fronts, and jumps can be of four types: shocks (or contact 
discontinuities), rarefaction waves, non-physical waves and zero- waves: 
J = 5U7^U7VPUZ. 

- Along each shock (or contact discontinuity) x^ = Xa{t), a € S, the val- 
ues of u~ = u'^''^{t,Xa—) and m+ = u'^''^{t, Xa+) are related by u"*" = 
''Pkai'^a)iu~) for some ka G {1, ...,n} and some wave-strength ctq. If the 
ka*^ family is genuinely nonlinear, then the Lax entropy admissibility 
condition < also holds. Moreover, one has 

I^Q - Afc„(u+,u")| < e 

where Afe^ (u"*", u~) is the speed of the shock front (or contact discontinuity) 
prescribed by the classical Rankine-Hugoniot conditions. 

- Along each rarefaction front — Xa{t), a €z TZ, one has u+ — (ctq)(u^), 
< da < e for some genuinely nonlinear family ka. Moreover, we have: 

\Xa - A/c„(m+)| < £• 
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- All non-physical fronts x ~ Xa{t), a 6 JVV travel at the same speed 
Xa — X > sup„ j |Ai(u)|. Their total strength remains uniformly small, 
namely: 

J2 \u''''{t,Xa+)-u'''\t,Xa-)\<e, yt>0. 

- The zero- waves are located at every point x — jh, j G ( — 7^, 7^) H Z. 
Along a zero-wave located at x^ = jah, a E Z, the values u~ — u'^'^{t, a;^— ) 
and = u'^'^{t, Xa+) satisfy — ^h{xon u~) for alH > except at the 
interaction points. 

- The total variation in space Tot.Var. •) is uniformly bounded for 
all t > 0. The total variation in time Tot.Var. x); [0, -|-cxj)} is 
uniformly bounded for x jh, j G Z. 

Finally, we require that ||u'^'''(0, .) — MoHlj^j-k-, — ^• 

Keeping h > fixed, we are about to let first e tend to zero. Hence we shall 
drop the superscript h for notational clarity. 

Theorem 3. Let v!^ he a family of e, h-approximate solutions of |2P^0^- There 
exists a subsequence u'^* converging as i ^ +00 in hj^^ ((0, -|-cxo) xM.) to a func- 
tion u which satisfies for any ip G C\ ((0, -l-oo) x M).- 

[uipt + f{u)ifx] dxdt 







/"OO / />/) 

/ Y.^{t,]h)[ / 



g [jh + s, u{t, jh-)] ds] dt = 0. (22) 



Moreover Tot.Var. u{t, •) is uniformly bounded and u satisfies the Lipschitz prop- 
erty 

\u{t',x)-u{t",x)\dx<C'\t' ~t"\, t',t">0; (23) 



Now we are in position to prove [T| Theorem 4] with our weaker hypotheses. 
As in [T] we can apply Helly's compactness theorem to get a subsequence u^' 
converging to some function u in Uioc whose total variation in space is uniformly 
bounded for all t > 0. Moreover, working as in [2| Proposition 5.1], one can 
prove that u'^'{t, •) converges in LI to u{t, •), for all t > 0. 

Theorem 4. Let u^^ be a subsequence of solutions of equation with uni- 
formly bounded total variation converging as i ^ -l-oo in LI to some function 
u. Then u is a weak solution to the Cauchy problem (Qp-flli- 

We omit the proofs of Theorem [3] and [3] since they are very similar to the 
proofs of Jj Theorem 3 and 4]. We only observe that, in those proofs, the 
computations which rely on the Loo bound on the source term have to be 
substituted by the following estimates. 
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Concerning the proof of Theorem [31 

rh 



\g{jh + s,u''{t,jh-)) - g{jh + s,u{t,jh~))\ ds 
< I \\g{3h + s,-)\\c^-\u'[t,jh-)-u{t,]h-)\ ds 



< Sh ■ \u%t,jh-) -uit,jh-)\ . 
• Concerning the proof of Theorem |3J 

\g{jh + s,u''{t,jh~))\ ds< Mjh + s, ■)\\ci ds < et 







and 



h 

h 



| g (j/i + s, u^{t, jh—)) — g {jh + s, u{t, jh + s)) | ds 

< / \\gijh + s,-)\\ci-\u''it,jh-)^u{t,jh + s)\ ds 
Jo 

< ih ■ Tot.Var. •), [{] - l)h, [j + l)h]] 



iu{x) lu'^t, x) — u{t, x) I . 

jh 

We observe that aU the computations done in [T, Section 4] rely on the source 
g only through the amplitude of the zero waves and on the interaction estimates. 
Therefore the following two theorems still hold in the more general setting. 

Theorem 5. There exists 6 > such that if 11(^11 ^^(K) sufficiently small, 
then for any (small) h > there exist a non empty closed domain VhiS) and 
a unique uniformly Lipschitz semigroup : [0, +cx3) x VhiS) — > T)h{5) whose 
trajectories u{t, .) = Pf^Uo solve and are obtained as limit of any sequence 
of e,h-approximate solutions as e tends to zero with fixed h. In particular the 
semigroups P^ satisfy for any Uo,Vo G Dhi^), t, s > 

P^uo = Uo, P^ o P^uo - PiVt^o: (24) 



for some L > 0, independent on h. 



+ t-s 



(25) 



Theorem 6. // ll^^ll^j^i-jf) is sufficiently small, there exist a constant L > 0, a 
non empty closed domain V of integrable functions with small total variation 
and a semigroup P : [0, +oo) x P ^ 2? with the following properties 



i) Pqu = m, Vu e V; Pt+sU ^ Pt o P^u, Vm G V, t, s > 0. 
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ii) \\PsU-Ptv\\j^^^^^<L(\s-t\ + \\u-v\\j^^^^^, VueV, t,s>0. 

Hi) for all Uq G T), the function u{t, •) = PtUo is a weak entropy solution of 
system flj). 

iv) for some 5 > and all h > small enough V C 'Dh{5). 

v) There exists a sequence of semigroups P^' such that Pj^'u converges in LI 
to PtU as i ^ +00 for any u E V. 

Remark 3. Looking at [U (4.6)] and the proof of [Ij Theorem 7] one reahzes 
that the invariant domains 'Dh{6) and 2? depend on the particular source term 
g{x,u). On the other hand estimate [U (4.4)] shows that all these domains 
contain all integrable functions with sufficiently small total variation. Since the 
bounds 0(1) in Lemma [5] depend only on / and on ||M||]^j^, also the constant 
Ci in [U (4.4)] depends only on / and on ||M||j^^. Therefore there exists S > 
depending only on / and on such that T>h{S) and V contain all integrable 

functions u{x) with Tot.Var. {u} < S. 



4 Uniqueness of BV entropy solutions 

The proof of uniqueness in [l] strongly depends on the boundedness of the 
source, therefore we have to consider it in a more careful way. 



4.1 Some preliminary results 

As in [1] we shall make use of the following technical lemmas whose proofs can 
be found in [3]. 

Lemma 6. Let (a, 6) a (possibly unbounded) open interval, and let X be an 
upper bound for all wave speeds. If u, v G T^hiS) then for all t > and h > 0, 
one has 

pb—Xt pb 

/ _ \{Ptu) (x) - (Pj^v) {x)\dx<L / \u{x) - v{x)\ dx. (26) 
Lemma 7. Given any interval Iq — [a,b], define the interval of determinacy 

It ^[a + \t,b- \t], t<^-^. (27) 

2A 

For every Lipschitz continuous map w : [0,T] ^ T>h{5) and h > 0; 

||z«(t)-PXO)||^^(,^^ (28) 

,r* r \\w{s + rj)-P;;wis)\\^ ] 
< LI niminf^^ ^ '^"^ }ds. 
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Remark 4. Lemmas [SI E] hold also substituting with the operator P. In 
this case we have obviously to substitute the domains "Dfi^S) with the domain 
V of Theorem [51 

Let now ug,Ur be two nearby states and A < A; we consider the function 

"<'.^)={;: ::>a;::, « 

Lemma 8. Call w{t,x) the self-similar solution given by the standard homoge- 
neous Riemann Solver with the Riemann data U5\) . 

(i) In the general case, one has 



+ 00 

\v{t,x) - W{t,x)\ dx = 0(1)|W£ - Ur\] (30) 



(ii) Assuming the additional relations Ur — Ri{a){ut) and A — Xi{ur) for some 
a > 0, i = I, . . . ,n one has the sharper estimate 

1 /■+°° 

-/ \v{t,x)-w{t,x)\dx = &{l)a'^; (31) 

^ J-oo 

(Hi) Let u* S and call A^ < . . . < A* the eigenvalues of the matrix A* = 
V/(m*). If for some i it holds A*(ur — ui) — A*(wj. — ui) and A = A* in 
I129\) then one has 



1 /•+°° / * * \ 

- / \v{t,x) — w{t,x)\dx — 0{l)\ui — Ur\i\ue — u*\ -\- \u* — Ur\); (32) 

^ J — oo ^ ' 

We now prove the next result which is directly related to our ft,-Riemann 
solver. 

Lemma 9. Call w{t, x) the self-similar solution given by the h-Riemann Solver 
in Xo with the Riemann data il5\) . 



(i) In the general case one has 

— j \v{t,x) — w{t,x)\dx — @{\)y\Ui — Ur\^ J uj{xo + s) dsj; (33) 

(ii) Assuming the additional relation 

-1 

Ur = Ui -\- [V/] (m*) / g{xo + s, u*)ds 
Jo 

with X — in i29\) one has the sharper estimate 
\v{t, x) — w(t, x) \ dx 

0(1) /" uj{xo + s) ds -\-\ue - u*\^ ■ I uj{xo + s)ds. (34) 



— oo 

h 
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Proof. Estimate is a direct consequence of Lemma \5\ Let us prove now 
H). Since A = we derive 



1 



\v{t, x) — w{t, x)\ dx 



1 



\uf — w{t, x)\ dx ~\ — / \ur — w(t, x)\ dx 



At 



At 



(35) 



0(1) 



.L — 1 L—p-\-l 



0(l)|fi| 



This leads to 



|6| = 



E [h, Ui, Ur] Xo] - E [h, ue, ^h{xo, uey,Xo] 

0(1) \Ur - <i>hiXo,Ue)\ . 



To estimate this last term, we define b{y, u) ^ f ^ ifi'^-) + v) and compute for 
some j/i, y2- 

ut + [Vfr\u*)yi-h{v2,ut) < 0(l)|2/i|-|"*-"f| + 0(l)|yi-2/2| 

+ [V/]""^ {ui)y2 -b{y2, ui) 

The function z{y2) — ui + [V/] ^ {ue)y2 — b {y2,ui) satisfies z(0) — 0, Dy2z{{)) 
0, hence we have the estimate 

ui + [V fy'^ {u*)yi-b{y2,uc) < 0(1) \yi\ ■ \u - ut\ + \yi - y2\ + \y2 

If in this last expression we substitute 

l*h I* h 

yi = / g{xo + s,u*)ds, y2 ^ g{xo + s,ui)ds 
Jo Jo 

then, we get 

\ur - ^h{xo,ue)\ ^ 0{1)(^ j uj{xo + s) ds + \ui - u*\J J Lj{xo + s)ds 
which proves ([M]) . □ 

4.2 Characterization of the trajectories of P 

In this section we are about to give necessary and sufficient conditions for a 
function u{t, •) £ P to coincide with a semigroup's trajectory. To this end, we 
prove the uniqueness of the semigroup P and the convergence of all the sequence 
of semigroups P'' towards P as ft, ^ 0. 



15 



We begin by introducing some notations: given a BV function u = u{x) and 
a point ^ G K, we denote by U^^,^-^ the solution of the homogeneous Riemann 
Problem ^T5\i with data 

ui = lim u{x), Ur — lim u{x), Xo = (36) 

Moreover we define U^^.^-^ as the solution of the linear hyperbolic Cauchy prob- 
lem with constant coefficients 

wt + Awx — g{x), w{0,x) = u{x), (37) 

with I = V/ (u(0), g{x) - g {x, u{0). 

We will need also the following approximations of U'^u-^) ■ f be a piecewise 
constant function. We will call w'^ the solution of the following Cauchy problem: 

/•'» 

iw'')t+A{w%=J2^i^~jh) g{jh + s)ds, w\0,x)^vix). 

Define m*=u(^) and let = Ai(u*), r.i = ri{u*), U — li{u*) be respectively 
the i*^ eigenvalue, the i*'' right/left eigenvectors of the matrix A. As in [T] w 
and have the following explicit representation 

w{t,x) = 'S2\{k,u{x - Xit)) + [ {li,g{x')) dx'\ ri 
w''{t,x) - f^|(/„tJ(a:-A,t)) + l(/,,G"(t,x))jr„ (38) 



where the function G^{t, x) = (t, x)ri is defined by 



E 

j: jh£{x — Xit,x) 

- E 

j: jhG{x,x — Xit) 



{li,g{jh + s)) ds if Ai > 



{k,g{jh + s))ds ifA, <0. 



Using ([3]) we can compute 

G':it,x)- 



{li,g{x')) dx 



x-Xit 



Hence, for any a, 5 G R with a < 6, we have the error estimate 



rb~\t 



I a+Xt 



\w{t,x) - w'\t, x)\dx< 0(1) 



\u{x) — v{x) \ dx + {b — a)eh 



(39) 



(40) 



(41) 



16 



From ((5^1) . it is easy to see that x) is piecewise constant with discon- 
tinuities occurring along finitely many lines on compact sets in the (t, x) plane 
for t > 0. Only finitely many wave front interactions occur in a compact set, 
and jumps can be of two types: contact discontinuities or zero waves. The zero 
waves are located at the points jTi, j G Z and satisfy 



W''{t,jh+) ~ W''{t,3h-) ^ [\/f]-' (U*) / g{jh + s)ds 



(42) 



Conversely a contact discontinuity of the i*'* family located at the point Xa{t) 
satisfies Xa{t) = Xi{u*) and 



w'''{t,Xa,it)+) - w''{t,xc,{t)-) = an{u*) 

for some ct £ R. 

Now, we can state the uniqeness result in our more general setting. 



(43) 



Theorem 7. Let P : V x [0, +oo) V he the semigroup of Theorem\^ and let 
A be an upper bound for all wave speeds. Then every trajectory u{t, •) — PiUq, 
Uq G T), satisfies the following conditions at every r > 0. 

(i) For every ^, one has 



lim — 

0~,o e 



u{r + e,x)-ul^^y^^^ {e,x) 



dx = 0. 



(44) 



(a) There exists a constant C such that, for every a < <b and < 6* < 
one has 



b-S\ 



a+e\ 



U{T + 0,X)-Ul^,y^^) {e,x) 



dx 



(45) 



< C 



-1 2 

Tot.Var. {u(r); (a, b)} + I oj{x) dx . 



Viceversa let u : [0,T] ^ T) be Lipschitz continuous as a map with values in 
il(R, M") and assume that the conditions (Oj), ImP hold at almost every time t. 
Then u(t, •) coincides with a trajectory of the semigroup P. 

Remark 5. The difference with respect to the result in p] is the presence of 
the integral in the right hand side of formula ()45|) . If lo is in Lcx), the integral 
can be bounded by 0(1) (5 — a) and we recover the estimates in [I]. Note also 
that the quantity 

^((a, 6)) = Tot.Var. {ii(r); (a, 6)} + f uj{x) dx 

J a 

is a uniformly bounded finite measure and this is what is needed for proving the 
sufficiency part of the above Theorem. 
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Proof. Part 1: Necessity Given a semigroup trajectory u{t, •) — PfU, u £ V 
we now show that the conditions {^j, {np hold for every r > 0. 

As in [T] we use the following notations. For fixed h, 9, e > we define 
Jt = J- u Jt° U Jf+ with 



J° — 



j: = 



e-(20-i + r)A,e-(i-r)A) ; 
e-(t-T)A,e+(t-T)A] ; 



(46) 



Let U^u(T)-i) piecewise constant function obtained from U^^^^^y^-^ {9, x) 

dividing the centered rarefaction waves in equal parts and replacing them by 
rarefaction fans containing wave fronts whose strength is less than e. Observe 
that: 

f + OC 



Applying estimate to the function U^^^^y^^ we obtain 



< L 



lim inf 



?7, 



f^, ,,Jt-T + n)- P!!UH ^ .,(t - t) 

(M(r);C)V i; (M(-r);5)V ^ 



Ll(,/t+„) 



(48) 



(it. 



The discontinuities of Uf'^, s ,^ do not cross the Dirac comb for almost all times 
t G {t,t + 9). Therefore we compute for such a time t: 

1 



1 



tuirm(t -r + V,x)- {Py^:^^y,^{t - r)) (x) 

Ufuirm^ -r + V,x)- {Pn'uli%);i)it - r)) (x) 



(49) 



dx. 



Define Wt the set of points in which U^^^^y^^{t — r) has a discontinuity while 
Zh is the set of points in which the zero waves are located. If rj is sufficiently 
small, the solutions of the Ricmann problems arising at the discontinuities of 

1 



yf\{t — t) do not interact, therefore 



U, 



(«(r) 



,^^^{t-r + ^,x)-(py^^^^y^^{t-T)) (x) 



dx 



E + E 

-r + V,y)- {PyiL\ry,,)it - r)) (y) 



dy 
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Note that the shock are solved exactly both in Uf''^, ,s and in P^Uf'^, . c\ there- 
fore they make no contribution in the summation. To estimate the approximate 
rarefactions we use the estimate (PT|) hence 



dx 



<0(1) E kP<0(l)eTot.Var.{c/;^(^)^^)(t-r);jo} (50) 



rare f action 



<0{l)e\u{T,^+)-u[T,^-)\ 

Concerning the zero waves, recall that t is chosen such that C/,"'^/ ^ ,^ is constant 
there, and is the exact solution of an /i-Riemann problem, hence we can 
apply (|33p with = Uj. and obtain 



Y. - 



Ulu(ry,i)it -r + il,x)- {Py^i^y^^^it - r)) (x) 



< 0(1) I ^{jh + s)ds<0il) I [ 



ijj{x) dx + Eh 



Finally using ((5T|) and ((50|) we get in the end 
1 



Ulu%y,oi*-r + ^,x) - (P,fC/f^.);«)(i-r)) (x) 



dx 



dx 



(51) 



(52) 



= 0(1) j y ^ tj(x) dx + ih + e j . 

Moreover, following the same steps as before and using ([50)1 and ([55]) with 

U£ = we get 



1 



V Jj+ 

t + ri 



UlU^)it-r + ^,x) - (p,;'C/«^^^^^)(t-r)) (x) 



(53) 



0(1) < y ^ w(x) dx + £/i 



Note that here there is no total variation of U^u{t)-C) ^i^*^^ ■^t it is constant. 
A similar estimate holds for the interval Jj^^. Putting together ([5^ . ([55]) . 
one has 



1 



'7 



<(.)^,)(t-r + ^,x) - [Py^^l^y^^^it - r)) (x) 

= 0(1)(/" c^(x) dx + e,, + e). 
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Hence, setting i = f/f;''(^).^)(0) = C/J„(^).^)(0) by (gH]), we have 



Finally we take the sequence P''* converging to P. Using we have 



< -L|mT)-C||Li(,,. 



(55) 



|u(t, x) — v{x) \ dx 



where eg tends to zero as 9 tends to zero due to the fact that u(r) has right and 
left limit at any point: for any given e > if is sufficiently small \u{t,x) — 
v{x)\ = \u{t,x) -uiT,^-)\ <eioTxe {^-2X9,^). 
Therefore by (gT]), jM]), we derive: 



1 



u{T + 9,x)-Ul^^y^^^ {9,x) 



Peu{T) - P,''-w(t) 



dx 



Li( 



+ eg + 0(1) 



uj{x) dx + ihi 



Jr 



The left hand side of the previous estimate does not depend on e and hi, hence 



^-ex 



u{t + 61, x) - C/»„(^).^,) (61, x) dx = 0(1) j uo{x) dx + eg. 

J T 



Note that the intervals Jr depend on 9 fsee So taking the limit as 6* ^ 
in the previous estimate yields 



To prove (EJ) let 6* > and a point (r, ^) be given together with an open 
interval (a, 6) containing ^. Fix e > and choose a piecewise constant function 
V ^ V satisfying v{£^) = w(t, ^) together with 




\v{x)-u{t,x)\ dx <e, Tot.Var. {i;;(a,6)} < Tot.Var.{u(r);(a,6)} (56) 



Let now w'^ be defined by §^ (u* = v{£,) = u(r,^)). From (gT]), ^ we have 
the estimate 

dx <0{l)(^e + eh{b-a)y (57) 



-8\ 
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Using (HIl), dMl) we get 



\w''{e,x)-{P^w''{0)){x)\dx (58) 

a+6\ 

rr+e Wwf^^t - T + ry) - P^'u,"(t - r)\\^,,, . 

< L / liminf-^^ ' ^ '+''' dt 

Jr ''^0 V 

where we have defined It+n = h-r+ir Let i G (r, t + 0) be a time for which 
there is no interaction in in particular, discontinuities which travel with a 
non-zero velocity do not cross the Dirac comb (this happens for almost all t). 
We observe that by the explicit formula ([38]): 



Tot.Var.|w''(t -T);/t| = 0(l)(Tot.Var.{w; (a,6)} + y oj{x) dx + ei}j (59) 

|w''(<-T,x)-t;(C)| =0(l)(Tot.Var.{w;(a,6)} + ^ uj{x)dx + eh^. (60) 
As before for -q sufficiently small we can split homogeneous and zero waves 

7'' (t - T + 77, X) - {Pl]w'' (t-T)) (X) \dx (61) 



1 



E + E \- \^\t-T + r^,x)-{P:;w\t~T)){x)\dx 



The homogeneous waves in satisfy (|43| . with w(^) in place of u*, hence 
we can apply (|32p which together with ((59|) . (|60p leads to 



-1 rX+ATj 

E -/ |w''(t-T + r;,x)- (P>''(i-T)) (a;)|dx 

< 0(1) ^ |Ai(;''(f-T,a;)|(Tot.Var.{w;(a,6)}+ / a;(a;)dx-|-s 

£ce/tnWt 

< 0(l)Tot.Var. |w''(i-r),/t| (xot.Var. {z;;(a,5)} + y uj{x) dx + £}}j 

< 0(l)(Tot.Var. {«;(«, 6)}+ f lj{x) dx + ef}j^ 



where /S.w^{t — t, x) denotes the jump of 'w'^{t — t) at x. 

The zero waves in satisfy ([42]) . hence we can apply (p4|) which together 
with (l60l) leads to 
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1 nX + ATj 

J2 - \w''{t^T + T],x)^{P^w''it^T)){x)\dx 

- X! / ujix + s) ds ■ (^Tot.Ya.T.{v;{a,b)} + J lu{x) dx + ih 



x£itr\Zu 



< 0{l) J uj{x) dx + ih^ (^Tot.Yar.{v; {a, b)} + J uj{x) dx + ih^ 

< 0(l)(Tot.Var.{u; (a,6)} + /" u}{x) dx + 



Let now P^' be the subsequence converging to P. Since w'^{Q) = v using 
(|57|) . ([55)1 . and the last estimates we get 



- u[T + e,x)-ul(,),^)[e,x) dx 

^ Ja+e\ 

\\P9u{t) - Pg'u{T)\\-^^^^^ ||M(r) - 

^ 7. h 1j 



+0(l)|£±£^ii_lL^+ (Tot.Var.{v;(a,6)} + ^'w(x) dx + ih^y 
So for e, hi — > we obtain the desired inequality. 

Part 2: Sufficiency By Remark|3]we can apply to P and hence the proof 
for the homogeneous case presented in [3] , which relies on the property recalled 
in Remark [5l can be followed exactly for our case, hence it will be not repeated 
here. □ 

□ 

Proof of Theorem [T] It is now a direct consequence of Theorems |6] and [71 □ 



5 Proofs related to Section [2] 

Consider the equation 

Ut + f{u)x = a'g{u) 

for some a G BV. Equation ^ is comprised in this setting with the substitution 
a I— > In a. For this kind of equations we consider the exact stationary solutions 
instead of approximated ones as in (|12p . Therfore call <i>(a, u) the solution of 
the following Cauchy problem: 

j^u{a) = [Duf{u{a))y^ 9{u{a)) f^r,^ 
u{Q) = u ^ ' 
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If a is sufficiently small, the map u i-^ $(a,u) satisfies Lemma [3] We call 
a-Riemann problem the Cauchy problem 



ut + f{u)x = a'g{u) 

{a-,ui) ifx<0 (63) 
(a+, Ur) if X > 



(a, w)(0, x) 



its solution will be the function described in Definition [T] using the map <i>(a+ — 
a~,u~) instead of the in there. Observe that if a+ — the a-Riemann 
solver coincides with the usual homogeneous Rieniann solver. 

Definition 3. Given a function u G BV and two states , a+, we define 
Ul{t,x) as the solution of the a-Riemann solver with u; — ""(0— ) and 

Ur — w(O-l-). 

Proof of Theorem [2j Since l|aj||j^j^ — |a+ — a^|, hypothesis (P2) is satisfied 
uniformly with respect to /, moreover the smallness of |a+ — a~\ ensures that 
the LI norm of w in (P3) is small. Therefore the hypotheses of Theorem [T] are 
satisfied uniformly with respect to 

Let be the semigroup related with the smooth section ai. By Remark 
[3l if Tot.Var. {u} is sufficiently small, u belongs to the domain of P' for every 
/ > 0. Since the total variation of P}u is uniformly bounded for a fixed initial 
data M, Helly's theorem guarantees that there is a converging subsequence P/'m. 
By a diagonal argument one can show that there is a converging subsequence of 
semigroups converging to a limit semigroup P defined on an invariant domain 
(see [U Proof of Theorem 7] ) . 

For the uniqueness we are left to prove the integral estimate p4)) in the 
origin with t/f subsituted by UK 

Therefore we have to show that the quantity 



^ p+ex 



-ex 



u{T + e,x)-Ul^^^{9,x) dx (64) 



converges to zero as 9 tends to zero. We will estimate (|64p in several steps. 
First define v = (r) ^) ^^^^ compute 

+ex 

\iPgu{T)){x) - {Pev){x) {9,x)\dx < ee. (65) 

-ex 

as in ([55|) . Then we consider the approximating sequence P'' corresponding to 
the source term a;, and the semigroups P'"'' which converge to P'' in the sense 
of Theorem |6l Hence we have 



lim lim — 

i^oo h^O 9 



+ex 

li.h 



-ex 



iPt''v)ix) - {Pev){x) 



dx = 
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For notational convenience we skip the subscript i in li. As in ([Tf|) we 
approximate rarefactions in C/'j.^-) introducing the function U^^^^y Then we 
define (see Figure [3]) 



U{x) 



for - ll2<x< 1/2 



Ul'^'^^(t~T,x-l) for x> 1/2 



where U{x) is piecewise constant with jumps in the points jh satisfying U{jh+) = 

jj\i,e,l,h 

X — Xt 




u{t, 0-) = ui 



u{t, 0+) = Ur 



X = I + At 



U 

Figure 3: lUustration of tj^'^^^'^ in the {t,x) plane 

^{jh, U(jh-)). Furthermore U{-l/2-) = Uf;^.^^{t - t,0-) and $ is defined as 
in p2p using the source term g{x,u) — a[{x)g{u). Observe that the jump be- 
tween U{l/2—) and t^f — t,1/2+) does not satisfy any jump condition, 
but as U {x) is an "Euler" approximation of the ordinary differential equation 
f{u)x — a'ig(u), this jump is of order Sh- Since U'^^^^ and C^fj'^'j''' have uniformly 
bounded total variation we have the estimate 



-ex 



ex 



U'ulr)i<^.x)-Ul;l;\\6,x) dx<0{l)- 
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the bound 0(1) not depending on h. We apply Lemnia[7]on the remainmg term 





+0\ 




ex 






< L 


i: 







dx 



To estimate this last term we proceed as before. Observe that P'^'' does not 
have zero waves outside the interval [— ^ — /i, ^ + /i] since outside the interval 
[— ^, ^] the function aj is identically zero. If 77 is small enough, the waves in 
^v'^^uIt)'^ ^ ''') '^'^ ^'-'^ interact, therefore the computation of the LI norm in 
the previous integral, as before can be splitted in a summation on the points 
in which there are zero waves in P ' or jumps in Uu(r)^{t-T). Observe that 

the jumps of C^f [t — t + rj) in the interval (— ^, + 5), are defined exactly as 
the zero waves in P'''' so we have no contribution to the summation from this 
interval. Outside the interval [— | — /i, | + /i], P^ coincides with the homoge- 
neous semigroup, hence we have only the second order contribution from the 
approximate rarefactions in U^^\^'^ {t — t) as in (|50p . Furthermore we might 

have a zero wave in the interval [— ^ — /i, — ^] and a discontinuity of U'^i^^^Y'' 
the point a; = ^ of order Sh- Using ([55]) for the zero wave and ([50)1 for the 
discontinuity (since P^ is equal to the homogeneous semigroup in a; ~ ^), we 
get 

liminf " "'"^ ^ " ^ < 0(1) + ) 

Which completes the proof if we let first e tend to zero, then h tend to zero, 
then I tend to zero and finally tend to zero. As in the previous proof, the 
sufficiency part can be obtained following the proof for the homogeneous case 
presented in [5]. □ 

Proof of Proposition [1] Call S the semigroup defined in [5] . The estimates 
for this semigroup outside the origin are equal to the ones for the Standard 
Riemann Semigroup see j3]. Concerning the origin we first observe that the 
choice ([9]) implies that the solution to the Riemann problem in [6, Proposition 
2.21 coincides with uK We need to show that 



I i-+e\ 
-0 ^ J -ex 



u{t + 9,x)- Ul^^^ {d, x) dx = 0. (66) 



with u{t,x) = {StUo){x). As before, we first approximate U^^^^ with U^^^^ 
and u{t) with Uf^(^^^{0)^v then we apply Lemma [7] (which holds also for the 
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semigroup S) and compute 



-e\ 



ex 



{Sev){x) 



IP'" 



{e,x) 



dx 



lim inf 



rj) 



[Jt- 



The discontinuities of U^'^^^ are solved by S*,, with exact shock or rarefaction 
for X ^ and with the a-Riemann solver in a; = therefore the only difference 
+ 77) and SjjUh^s {t — r) are the rarefactions solved in an 



between t/'f > (t 



approximate way in the first function and in an exact way in the second. Recall- 
ing pip we know that this error is of second order in the size of the rarefactions. 
To show that ([66|) holds, proceed as in (|50| . □ 
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